Strange hadron resonances as a signature of freeze-out dynamics 
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We study the production and the observability of A* (1520), 7f* (892), and £* (1385), strange hadron 
resonances as function of the freeze-out conditions within the statistical model of hadron production. 
We obtain an estimate of how many decay products are rescattered in evolution towards thermal 
freeze-out following chemical freeze-out, and find that the resonance decay signal is strong enough to 
be detected. We show how a combined analysis of at least two resonances can be used to understand 
the chemical freeze-out temperature, and the time between chemical and thermal freeze-outs. 
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Hadronic particle signatures of the formation of QGP 
phase in relativistic heavy ion collisions are most sensitive 
when the dense hadron matter fireball breakup is sudden 
But final state particles could also emerge remem- 
bering relatively little about their primordial source, hav- 
ing been subject to rescattering in purely hadronic gas 
phase |jP]. Which reaction picture applies can have de- 
cisive influence on our understanding of the underlying 
physics, and thus we propose here a systematic method 
to experimentally make a distinction. 

To address this question we consider strange hadron 
resonance abundances. At this time A(1520) has been 
observed in heavy ion reactions at SPS energies fol- 
lowing a suggestion that such a measurement was possi- 
ble @]. SPS H and RHIC experiments §] report mea- 
surement of the if* (892) signal, and RHIC has already 
measured both the K* a and the K*°. In the SPS case, 
the A(1520) abundance yield is about 2.5 times smaller 
than expectations based on the yield extrapolated from 
nucleon-nucleon reactions. This is to be compared with 
the enhancement by factor 2.5 of A-production in the 
same reaction in terms of the same comparison. 

A possible explanation for this effective suppression 
by a factor 5, or more, is that the decay products (tt, A) 
produced at a high chemical freeze-out temperature T ~ 
175 MeV have rescattered and thus their momenta did 
not allow to reconstruct this state in an invariant mass 
analysis. However, the observation of a strong if *°-yield 
signal contradicts this point of view, since the if* (892) 
decays faster (r K ,o = 50.8 ±0.9 MeV > r A . = 15.6 ± 1 
MeV). Therefore, if* (892) should be even more sup- 
pressed: a back of envelope calculation based on expo- 
nential population attenuation suggests that if the ob- 
servable yield of A(1520) is reduced by factor 5, the ob- 
servable yield of if *°(892) should be suppressed by a fac- 
tor 10. 

Another explanation is that the chemical freeze-out 
temperature which governs the production of these res- 
onances in a thermal model is considerably lower, and 
hence rescattering of decay products is less significant. 
This is a point of view arising from a recent analysis of 
the hadronization process jjj . 

We explore the production, and the suppression of the 



observability of these resonances(if *°(892), A(1520)), 
and also explore the (more difficult) measurement of the 
E*(1385) state as a further way of distinguishing between 
reaction scenarios. We will show in a quantitative analy- 
sis how these measurements can constrain both the chem- 
ical freeze-out temperature and the lifetime of the hadron 
interacting phase evolving between the chemical and the 
thermal freeze-out conditions. 

The E*(1385) is expected to be produced in a thermal 
model more abundantly than A* (1520) in a hadronic fire- 
ball due to it's high degeneracy factor and smaller mass. 
Because of its 3 times shorter lifetime (Ts* = 36-39 
MeV > r A « = 15.6 ± 1 MeV), the E* signal is more 
strongly influenced by final state interactions than that 
of A*(1520). Like for K*°(892), one would naively expect 
that the observable yield of E* should be suppressed by 
a factor 10. 

We will first show that, in a fireball hadronizing ac- 
cording to the statistical model, the proportion of A s 
and Kaons resulting from decays of the main resonances 
can be noticeably dependent on chemical freeze-out tem- 
perature. Hadrons produced directly from a medium at 
temperature (T) much lower than the particle mass (as 
in all cases considered here) fill the available statistical 
phase space which has the relativistic Boltzmann distri- 
bution shape: 
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Here g is the statistical degeneracy, and 7; are 
the fugacity and equilibrium parameters of each valence 
quarks, and E is the energy. When the fireball is ex- 
panding at a relativistic speed, equation [I] describes the 
energy distribution of an element of the fireball in a ref- 
erence frame at rest with respect to the expansion (flow) . 
However, in this paper we evaluate ratios of particles with 
similar masses, and interaction modes, often considered 
in full phase space. For this reasons to a good approxima- 
tion, the flow effects largely cancel out. Similarly, for ra- 
tios of particles with the same valence quark composition, 
such as E*/A, A(1520) /A and, in the limit of X u = X d (as- 
sumed hence forward), if *° (892)/Jf " (= K*°/K + ) the 



1 



chemical factors (A's and 7's) cancel out between the two 
states compared. 

We now include in the yields the descendents arising 
from resonance decays. In a decay of the type R — > 1 + 2 
where a particle R (mass M, transverse mass Mt, rapid- 
ity Y) into a particle 1 (mass m, transverse mass my, ra- 
pidity y) and 2 (mass m2), the distribution d 2 N\ / dm\,dy 
can be obtained in terms of d 2 Nn/dM^dY using the fol- 
lowing transformation : 



dN x 
dm^pdy 

J = 



g r b 

Airp* 



C f MT+ 2 d 2 N R 

/ dY / dMlJ —3 " 
y- Jm t - dM 2 dY 



IM T - 
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Here AY = Y — y, y/s is the combined invariant mass of 
all the decay products except 1, and E* = ^jg(M 2 — to 2 — 
m 2)> P* = V E* 2 — m 2 are the energy and momentum of 
decay particle 1 in the resonances rest frame. g r accounts 
for a resonances degeneracy, while b is the branching ratio 
of the considered decay channel. 

A full justification for Eq. (§) is found in Refs. gQ, 
where it is obtained by transforming the resonance dis- 
tribution from it's rest frame into the fireballs rest frame, 
and integrating over all the kinematically allowed values 
of rapidity and transverse mass. J is the Jacobian of this 
transformation, with integration limits: 

Y± = y±smh _1 (— ) 
rriT 



In Fig. [l] we show the relative thermal production ra- 
tios at chemical freeze-out over the entire spectrum of 
rapidity and tot (solid lines) as well a central rapid- 
ity range defined by the y — mj region covered by the 
WA97 experiment (\y\ < 0.5 in the center of mass frame) 
(dashed lines). The opportunity to measure the chemical 
freeze-out temperature by a measurement of the relative 
resonance yields is apparent. 

For example, at the lowest current estimates (T ~ 100 
MeV) of the final break up temperature in 158A GeV 
SPS collisions 33% of A s are actually primary E* s, the 
percentage rises to slightly more than 50% if chemical 
freeze-out occurs at T = 190 MeV. We recall that the 
CERN WA97 collaboration has measured hyperon yields 
to considerably greater precision than required to distin- 
guish between these limiting cases (jn], and the same is 
certainly expected of RHIC experiments such as STAR. 

Considering the practical feasibility of such a measure- 
ment it should be noted that the decay E* — > Ait would 
appear as a A-7T pair arriving from the central fireball. 
Experiments usually reconstruct S — > Att decays by find- 
ing the invariant mass of Att pairs with a common origin 
outside the collision, so a E* event would normally be 
discarded as a pair of unrelated particles coming directly 
from the fireball. An invariant mass analysis with res- 
olution better than E*(1385) natural width (35MMeV) 
would be needed to distinguish it from from S(1321) and 
minimize the combinatoric background. 
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E*m T cosh(AY) ± p T \ P* 2 - m 2 sinh 2 (AY) 

M± = M ^ 

simi (AY") 

Table | summarizes the dominant decay processes con- 
sidered in our analysis and their parameters (Clebsh- 
Gordon coefficients have been used to estimate decays 
such as (N*° -> A+7r-)/(A*° -> A tt )). 



TABLE I. A production mechanisms 



g 


Reaction 


P* 
MeV 


branching 
ratio 


visible? 


ps 4 


E* u (1385) -> EV J 


127 


^4% 


No 


8 


E* ± (1385) -> Att* 


208 


88% 


Yes 


4 


£* u (1385) -> Att u 


208 


88% 


No 


2 


S u - A 7 


74 


100% 


No 


4 


A(1520) -> NK 


244 


45% 


Yes 


3 


K* u (892) -> K+tt- 


291 


67% 


Yes 


3 


^i(1270) Kp 


76 


41% 


No 


3 


ifi(1270) -> K*"(892)n 


301 


16% 


No 


6 


i^* u (1400) -» A"* u (892)tt 


301 


100% 


No 


6 


_K"* U (1400) A"* u (892)tt 


301 


100% 


No 


5 


K*"(U30) -> A"* u (892)7r,7T7r 


« 400 


40% 


No 


5 


^* u (1430) -» Kn 


622 


50% 


No 


1 


K (1430) -> Ktt 


622 


100% 


No 
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FIG. 1. Temperature dependence of ratios of E*, K*° and 
A(1520) to the total number of observed K + A s,S s and Q s. 
Branching ratios are included. Dashed lines show the result 
for a measurement at central rapidity Ay = ±0.5. 

The momentum of each particle in the E* — > Att 
pairs CoM (center of momentum) frame would be « 208 
MeV/c, very close to the 211 MeV/c of the Q -> AK de- 
cay. As Fig. [I] shows (top two solid/dashed lines), the po- 
tentially observable (i.e. produced) E* s should be more 
abundant than these two hyperons ( 5, CI) at all consid- 
ered hadronization temperatures. The chemical poten- 
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tials required to compute the two ratios E*/f2,E*/S in 
Fig. [I] were taken from J3] . 

The ratios of observed particles, however, can be con- 
siderably different from production ratios, since if the 
decay products rescatter before detection their identifica- 
tion by reconstructing their invariant mass will generally 
not be possible. While the lifetime of the S and f2 are 
large enough to ensure that only a negligible proportion 
of particles decay near enough to the fireball for rescatter- 
ing to be a possibility, the lifetime of AT*°,E* and even 
A(1520) is within the same order of magnitude of the 
fireballs dimensions i.e. 2R/c ~ f/r. For this reason, 
a considerable number of decay products will undergo 
rescattering, and the estimation of this percentage is re- 
quired before any meaningful parameters are extracted 
from the data. 

Note in Fig.jl] that the relative E*/2 signal is remark- 
ably independent (within 5 %) of Temperature. This 
feature arises because the S*(1530) contribution cancels 
nearly exactly the thermal suppression of the S originat- 
ing in the 5 — E* mass difference. This effect could be 
used for a direct estimate of the E* lost through rescat- 
tering, even without knowledge of the freeze-out tem- 
perature, should the chemical parameters (A^ and 7,) be 
known. A simple test of hadronization model consists 
in measurement of the ratio E*/S. If it is significantly 
smaller than unity, we should expect a re-equilibration 
mechanism to be present. Otherwise sudden hadroniza- 
tion applies. 

We can go further in the use of the suppression of 
the considered resonances as a tool capable of estimating 
conditions at particle freeze-out. To do this we need to 
estimate the effect rescattering would have on the res- 
onances signal. We formulate a simple model based on 
the width of the resonances in question and the decay 
products reactions within an expanding fireball of nuclear 
matter. While more realistic models, possibly involving 
event generators are needed for precise quantitative anal- 
ysis, the model described below allows us to make a qual- 
itative prediction and to verify the feasibility of the S* 
measurement. 

We consider the decay of a generic resonance N* 



N* 



P1P2 , 



(3) 



in a gas of pions and nucleons. For the invariant mass 
of N* to be undetectable, it is assumed here that it is 
sufficient that either of the generic decay products p\ or 
P2 undergo one interaction. This interaction depends on 
the cross-section of the decay product with each parti- 
cle in the fireball (pions, nucleons and antinucleons), the 
speed of each decay product relative to a typical fireball 
particle and the pion and nucleon density in the fireball 
(decreasing with time for a fireball constantly expanding 
with flow v). 

This reaction rate is 



Pi = (cIttPtt + oinRn + ^iivPlvX 



I? 



R + vt' 



(4) 



and analogously for the second decay product '2'. Here, 
ct's are the (energy averaged) interaction cross sections, 
the p's are the densities, R is the fireball radius at chem- 
ical hadronization, (v) the velocity factor and v is the 
space-average of the flow velocity. 

To a good approximation, the thermal velocities in the 
fireball average out and (vj)i.2 = p*/ m i-2 ; here p* is as 
defined in Eq. (|^). In addition to scattering and decay 
we also in principle can allow that E* and it each have a 
probability of escaping the fireball. We have studied this 
effect and in this study of the measureability of the life- 
time of the interacting phase, such an escape probability 
is not relevant. The population equation describing the 
scattering loss abundance (N,) are, therefore: 



dNj 
dt 



-N N > - NiPi 



1,2 



dNff* 1 

— ^— _ --N N , 



(5) 



The required nucleon density pat is obtained in the a 
relativistic Boltzmann approximation: 



Pn 



fJ 



A^{\ qlq fTK 2 (-) 



(6) 



where the K2(x) is the 2nd modified Bessel function. 
We consider the nucleons to have a mass of ~1 GeV, 
and a degeneracy of 6, to take the p,n and the ther- 
mally suppressed but higher degeneracy A contributions 
into account in a first approximation. Considering the 
CERN/SPS environment we calibrate the fireball size R 
to obtain the baryon multiplicity to be 360 within the 
volume, and we take R = 8 fm-145/T[MeV]. 

The temperature scaling of R is chosen to assure that 
with the pion density computed in the massless particle 
limit, leading to 



An 2 
~90 



(7) 



the entropy of the expanding system is conserved. The 
effective degeneracy A is chosen to fix the entropy per 
baryon to be equal to the observed number of ~40. We 
believe that this qualitative model of baryon and meson 
density has the right physical magnitude and scales cor- 
rectly with variation of freeze-out temperature. 

To estimate the average cross sections needed in 
Eq. (ji|) we recalled that the NA49 collaboration using 
the UrQMD event generator pSfl , has found that « 50% 
of A(1520) undergo rescattering after a time of 80 fm/c. 
In order to reproduce this result at T « 160 MeV, we 
took scattering cross sections shown in the top portion 
of table ||, based on Refs. [p^|Jl4[ . It turns out that the 
cross section on nucleons and antinucleons are sufficient 
to make the fireball essentially opaque to one of the res- 
onance decay products. 
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FIG. 2. Produced (dashed line and observable (solid lines) 
ratios A(1520)/(total A) (top) and S* /(total A) (middle) and 
K /K (bottom). The solid lines correspond to evolution 
after chemical freeze-out of 1,2,3,4,5,7,10,15,20 fm/c, respec- 
tively. 

The uncertain meson scattering cross sections (denoted 
in Table |l| for irir by X and for ttK by 0.5 • X) are there- 
fore not material for the results we present. We quan- 



titatively explored the range 1 < X < 100 mb without 
seeing any significant change in our results. In the p res- 
onance region applicable at the conditions considered the 
value X ~ 70 mb is appropriate H] . Remainder of Table 
presents an overview of other parameters used in the 
model. 

Figure|| shows the dependence of the A(1520) /A, E*/A 
and K* (892) / K + on the temperature and lifetime of the 
interacting phase. It is clear that, given a determination 
of the respective signals to a reasonable precision, a quali- 
tative distinction between the high temperature chemical 
freeze-out scenario followed by a rescattering phase and 
the low temperature sudden hadronization scenario can 
be made. We also note that despite the shorter lifetime 
of the E* and higher pion interaction cross section, more 
E* decay products should be reconstructible than in the 
A* (1520) case, at all but the highest temperatures under 
consideration. 

Fig. H demonstrates the sensitivity of strange hadron 
resonance production to the interaction period in the 
hadron phase, i.e. the phase between the thermal and 
chemical freeze-out. What we learn from Fig. ^ is that 
while the suppression of one of these ratios considered 
has generally two interpretations, as it can mean either a 
low temperature chemical freeze-out or a long interacting 
phase with substantial rescattering, the comparison of 
two resonances with considerably different lifetimes can 
be used to constrain both the temperature of chemical 
freeze-out and the lifetime of the interacting phase. 

We can reexpress the results presented in Fig. || repre- 
senting one ratio against the other as is seen in Fig. and 
attempt a first comparison with experiment. Combining 
the K*°/tt+ results presented in with earlier K~ Jit" 
and K + /n + results £j]|, it seems that the K*° / K + ra- 
tio should be w 0.2. Unless the A(1520)/A ratio will 
be considerably different than what is expected from [gj , 
experimental measurement for SPS energies should be in 
the top left corner of the bottom portion of the Fig. ||, 
corresponding to a low chemical freeze-out temperature 
(w 145 MeV) and a short, if not negligible, time spend 
in the hadron re-interaction phase. Pending more exper- 
imental work, this evaluation serves only to demonstrate 
the functioning of the method we propose. 



TABLE II. Scattering model parameters 



ovjv (mb) 


0~KN 




0VX 


CNN 


a NN 


24 


20 


X 


0.5- X 


24 


50 


r_. 


Ta(1520) 


Pit *°(892) 


35 MeV 


15.6 MeV 


50 MeV 


escape rate (fm 6 ) 


negligible 


v 


0.5 


R(fm) 


8145/T [MeV] 


Mb 


220 MeV 


A 


11 
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Should the E*/S ratio become available along with a 
determination of Xi and ji, both the temperature and the 
lifetime can be inferred from the E* alone. Fig. |J shows 
an analogous diagram to the one shown previously, but 
with the E*/(all A) on the horizontal axis and E*/(all 
H) on the vertical. Since E*/(S + S*) only depends on 
temperature by about 5%, see Fig. |l], the percentage of 
rescattered E* can be, to a good approximation, read off 
the vertical axis in Fig.H. 
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FIG. 3. Dependence of the combined A(1520)/(all A), 
E*/(all A) and K*° (892) /(all K-) signals on the chemical 
freeze-out temperature and interacting phase lifetime. 

We have shown that a comparison of two or better 
three observable strange hadron resonance abundances, 
with different masses, widths and interaction modes of 
the decay products, can be used to estimate the cur- 
rently uncertain hadron freeze-out behavior. We have 
chosen to consider strange hadrons as they have rela- 
tively narrow widths. A similar study of the lifespan 
between thermal and chemical freeze-outs, and chemical 
freeze-out temperature as parameters, can be undertaken 



with non-strange nucleon resonances. We believe that 
the experimental difficulties of such a measurement are 
much greater: the relative yield of resonances such as 
AT* (1440), A(1230) compared to the nucleon is smaller 
than in case of strange baryons, the widths are much 
larger (Fat. ~ 300 MeV, T A = 120 MeV), and the final 
state nucleon has a greater scattering cross section in 
hadronic matter than the strangeness carrying hadron. 
Therefore the detectability of the non-strange baryon res- 
onances in terms of an invariant mass analysis in the de- 
cay channels should be much smaller, except if a (very) 
sudden hadronization applies. Seen from that perspec- 
tive, it will be interesting to see if the A(1230) state can 
be observed at all, as this would be suggesting that chem- 
ical and thermal freeze-out are indeed nearly coincident. 
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FIG. 4. Dependence of the combined E*/(all A s) and 
E*/(all E s) signals on the chemical freeze-out temperature 
and interacting phase lifetime 

We believe that a detector capable of reconstructing 
5 and f2 decays, with a narrower invariant mass resolu- 
tion than the natural E*-width (35 MeV) may be capable 
of measuring the E* signal by performing an invariant 
mass fit on A-7T pairs arriving directly from the fireball. 
Present day state of the art detectors such as STAR and 
NA57 should satisfy this criterion. We have shown that 
the E*/S ratio is remarkably independent (within 5 %) of 
freeze-out temperature and is greater than unity. A sig- 
nificantly depleted ratio would directly measure rescat- 
tcring of hadrons after chemical freeze-out. Considering 
other ratios we learn from Fig. [2] that while the suppres- 
sion of one of particle ratios considered is generally not 
able to allow resolution of the freeze-out properties, as it 
can mean either a low temperature chemical freeze-out 
or a long interacting phase with substantial rescatter- 
ing, the comparison of two relative resonance yields with 
considerably different lifetimes can be used to constrain 
both the temperature of chemical freeze-out and the life- 
time of the interacting phase. This is than qualitatively 
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evaluated in Fig. |[ 

We have shown in quantitative fashion how a measure- 
ment of the relative yields of A(1520) and A'*°(892) with 
respect to the ground states A and K allows to constrain 
the chemical freeze-out temperature, and the lifetime of 
the intermediate evolution period ranging between ther- 
mal and chemical freeze-outs in relativistic heavy ion 
collisions, and hence to distinguish between sudden and 
staged hadronization scenarios. We have suggested that 
present experiments should be also able to detect the 
S* resonance, which when subject to the same analysis 
method offers a consistency check for the reaction mech- 
anism considered. 
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